Magnetic field relaxation in ferromagnetic Ising systems 



R. Melin 

CRTBT-CNRS, 38042 Grenoble BP 166X cedex France 

Abstract 

We analyze the thermal magnetization reversal processes in magnetic grains. 
Two experiments are carried out: swtiching time and switching field experiments. 
In both cases, we find that the simulated behavior is coherent with existing exper- 
imental data (the streched exponent of the switching time experiment increases 
with the temperature and is superior to unity; there exists a master curve for the 
switching field experiment). Moreover, we simulated magnetic grains in a region of 
parameters where no experimental data are available. We find that the relaxation 
time distribution P(lnr) is gaussian, and we find the existence of a strong field 
regime. 
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Magnetic grains are considered as good candidates to bring to the fore macroscopic quan- 
tum tunelling |l] . Before considering the problem of quantum tuneling at low temperatures, it 
would first be interesting to understand thermally activated processess at higher temperature. 
The aim of the present paper is to address this question in a simplified model. Experimentally, 
small ferromagnetic grains were studied for many years in powders, so that data were avail- 
able only for a collection of interacting magnetic grains, with different sizes and compositions. 
Using magnetic force microscopy and Lorentz microscopy ||, it was possible to resolve 
the magnetic properties of individual single-domain particles. The present article is based on 
the experimental works [|J, where it was possible to study the magnetization reversal of a 
single monodomain magnetic particule of size 0.1/zm x 0.05/im. Due to dipolar interactions, 
anisotropy aligns the spins along the largest direction of the grain. The spins are essentially 
XY spins since the grain is essentially bidimensional. We do not carry out direct simulations 
of XY spins at finite temperature, but, instead, we consider a simplified Ising spin model on a 
disk to simulate a magnetic grain. This means that a direct comparison with the experimental 
exchange couplings and magnetic field is not possible. However, we believe that our simula- 
tions are qualitatively valid, eventhough they do not directly represent magnetic grains. The 
reason why we choose the disk geometry is that the frontier is expected to play a great role 
in the magnetization reversal processess. In || , the boundary conditions are periodic, so that 
boundary effects are not properly taken into account. In order to compare with existing ex- 
periments |4[], we use the same experimental techniques than in [[|. Two types of experiments 
are possible 1) switching time experiments 2) switching field experiments. 

In the first class of simulations, the system starts in a low temperature configuration 
where nearly all the spins are paralell. A constant magnetic field opposite to the initial 
magnetization is applied from time t — 0. The initial configuration, antiparalell to the field, 
is then metastable, and the magnetization reverses to the oppposite magnetization. The time 
r corresponding to the cancellation of the magnetization is recorded. Since the reversal of 
magnetization process is stochastic, r is a random variable. In our simulations, we took 



1 



into account 30000 different runs. This number is necessarily huge in order to have a good 
statistics. In existing experiments, data with such an important number of reversal processes 
are not available at the moment, so that some of the results that are reported here are not 
directly comparable to experimental data because of the lack of experimental data in these 
temperature and field ranges. We first show that the distribution P(lnr) of the logarithm of 
the relaxation times is gaussian. This result is in agreement with other theories, especially 
for spin glasses |J, where a gaussian distribution of the logarithm relaxation times has been 
assumed, in the context of a phenomenological renormalization. We also analyze the mean 
value (lnr), and the variance ((lnr— (lnr)) 2 ) as a funtion of temperature. We find that, if the 
field is not too strong, both the mean value and the variance diverge as T approaches zero. If 
the field is too strong, the divergence is replaced by a saturation of these two quantities, so that 
long relaxation times at low temperature are inhibited by a too strong magnetic field. This 
qualitative behavior between a strong magnetic field regime and a weaker magnetic field regime 
is to be compared with the results of ||, where a strong magnetic field regime is exhibited. 
In the switching time experiments, we can study the variations of (In r) as a function of the 
variance ((lnr — (lnr)) 2 ). Within the context of phenomenological spin glass theories 0, this 
curve is expected to be linear. However, we do not find linear variations. Depending on the 
strength of the magnetic field, the curvature is negative (strong magnetic fields) or negative 
(weak magnetic fields). 

The second possible class of experiments is the so-called switching filed experiments. The 
idea is to start from a low temperature configuration, and to switch a field oppposite to the 
global magnetization. The field increases linearly with a rate v = —dH/dt. In the presence 
of the external field, the magnetization reverses and one measures the magnetic field H sw for 
which the magnetization cancels. The switching field H sw is averaged over a large number 
of runs, and computed for various temperatures. In order to analyze the data, we use the 
phenomenological model of reference 0. We find that, as in real experiments, there exists a 
master curve for (H SW )(T, v). 

1 Switching time experiment 

1.1 Square lattice with free boundary conditions 

In order to emphasize the importance of the boundary conditions, we first consider the case 
of an Ising model on a square. In this case, the domains of reversed spins can nucleate 
from the corner. The barrier to nucleate a domain of spins paralell to the magnetic field is 
Ei(R) = ttR(J — hR/2), where we have assumed that the domain of reversed spins is circular 
(this is an approximation, since the droplets are not necesserily circular). The maximum over 
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R of E\ (R) is nothing but the barrier that need to be paid to reverse the magnetization. This 
barrier is Ei(R[) = -nJ 2 j2h. In order to compare with the absence of edges (Ising model on a 
disk), we calculate the barrier to nucleate a domain of spins paralell to the field on the disk. 
We assume that the radius of the disk is large compared to the radius of curvature of the 
droplet. The barrier is given by -E^-R) = 2ttR(2J — hR). The first term AirRJ is the energy 
term associated to the antiparallel links at the boundary and —2irhR 2 is the magnetic energy. 
E2(R) is maximum for i?2 = J/K which corresponds to the critical size of the domain. If 
R < R2, the domain has a trend to regress whereas it naturally grows if R > R^- The barrier 
height is ^(-R^) = 27rJ 2 //i. We see that the barrier in the presence of edges is one half of 
the barrier in the absence of edges, from what we conclude that the spin system on a square 
with open boundary conditions is very sensitive to the presence of edges and will first flip 
from the edges, at least if the magnetic field is not too strong. To check this assertion, we 
took a snapshot of the system during a given run. The result is pictured on figure [I]. It is 
clear that the edges play an important role in the magnetization reversal, since the droplets 
are pinned at the edges, where the energy barrier is smaller than in the bulk. The pinning 
of the droplets is here purely geometric. One could also imagine that the droplets are pinned 
due to the presence of impurities. 

In what follows, we aim to elimitate the corner nucleation of domains. This is why we 
consider an Ising model on a disk, with free boundary conditions. 

1.2 Distribution P(lnr) 

We call r the time corresponding to the cancelation of the magnetization. We plotted on figure 
^| the distribution P(lnr) as a function of lnr. We see that these distribution in very good 
agreement with gaussian distributions. This means that the distribution of relaxation times 
in an external magnetic field is gaussian, provided one uses the logarithm of the relaxation 
times. This result is in agreement with the assumption of ||. In the context of spin glasses, 
it is assumed that P(lnr) is Gaussian. However, as we will see later, not all the assumptions 
of this phenomenological theory are valid in our model. 

1.3 Mean value and variance of the relaxation time distribution 

Since P(lnr) is a well defined gaussian distribution, we can study the variations of (lnr) and 
of the variance ((lnr — (lnr)) 2 ) as a function of temperature. We first begin with the mean 
value (In r) . The results for three different magnetic fields are plotted on figure |3|. We observe 
that if the intensity of the magnetic field decreases, the relaxation times increase. Moreover, 
if there exists a critical value h* f the magnetic field such that if \h\ < h*, the relaxation times 
diverges in the zero temperature limit, whereas if \h\ < h*, the relaxation times saturate in the 



3 



zero temperature limit. We identify this behavior with the existence of the strong magnetic 



The variations of the variance ((lnr — (lnr)) 2 ) as a function of temperature are plotted on 
figure |]. These variations are qualitatively similar to the variations of the mean value: if the 
intensity of the magnetic field is too strong, we observe a saturation of the variance, whereas 
it diverges for weaker magnetic fields. 

We are now in position to study the variations of the mean value (In r) as a function of the 
variance ((lnr — (lnr)) 2 ). These variations are pictured on figure [|. We observe that in the 
strong magnetic field region, the variance and the mean value are quite small, and both tend 
to a constant if the temperature decreases to zero. By contrast, we observe that if h = — 1.5, 
the curvature is negative, whereas it is positive if h = — 1. It would be interesting to know 
whether experimental on magnetic grains could confirm these simulations. However, these 
data are not available at the moment, due essentially to the fact that it is necessery to repeat 
the experience a large number of times to obtain a good statistics, and also to the fact that 
the experiments are not carried out for a given magnetic field, in order to avoid very long 
relaxation times. If we compare our results to the prediction of spin glass theory ||, we do 
not find that (lnr) is linear as a function of ((lnr — (lnr)) 2 ), and, according to the strength 
of the magnetic field, we find a positive or a negative curvature. We conclude that our model 
for thermal magnetic grains shares some common features with spin glasses since P(lnr) is 
found to be gaussian in both cases. However, whereas spin glasses phenomenological models 
predict a linear variation of (lnr) as a function of the variance ((lnr — (lnr)) 2 ), we do not 
find here such a linear behavior. 

1.4 Streched exponent (3 

An other possible analysis of the data of the switching time experiment if to buid the histogram 
P(r), and to fit the integrated histogram by a streched exponential 



This analysis is carried out on experimental data in [|J, where the streched exponent (3 is 
studied as a function of temperature. The main results are that (3 increases with tempera- 
ture, and may be superior to unity. As we shall see, our simulations agree with these two 
experimental facts. At this point, it should be underlined that the streched exponent j3 has 
nothing to do with the streched exponent of || in the context of spin glasses. In [|J, the 
autocorrelation functions of spin glasses are fitted by a streched exponential, and a streched 
exponent is deduced from the numerical data. In this case, the streched exponent is necessery 
inferior to unity, whereas in what we study, the streched exponent may be larger than unity. 



field regime of reference ||. 




(1) 
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We show on figure |^ the variation of the streched exponent (3 as a function of T, for two 
magnetic fields. It is clear that (3 increases with T, and may be superior to unity. We also 
observe that the streched exponent increases if the strength of the magnetic field increases, 
which is coherent with the fact both the mean value (lnr) and the variance ((lnr — (lnr)) 2 ) 
decrease if the strengh of the magnetic field increases. 

2 Switching field experiment 

In the switching field experiment, a linearly increasing magnetic field is applied opposite to 
the direction of the initial magnetization, with a rate v = —dH/dt. The magnetic field H sw 
corresponding to the cancellation of the magnetization is measured. In order to analyze the 
data, we follow reference [ffl, and use the phenomenological model of M. In this model, the 
average switching field (H sw ) is given by 

(H sw ) = H° sw (l-f a (T,v)), (2) 

with 

Hg W , Eq and c are some constants, a - 

e = {l-H sw /Hiy/«. (4) 

The experimental data consist of (H sw )(T,v) as a function of the temperature T and the 
switching rate v. The relation (0) predicts the existence of a master curve for (H sw )(T,u). 
We carried out numerical simulations for the square lattice on the disk, with open boundary 
conditions. We come to the conclusion that there exists a master curve for (H sw )(T,i/). 
The master curve is plotted on figure [7]. We see that there exists differences between the 
master curve and the behavior (H sw )/H® w = 1 — f a (T, v). These deviations also exist in the 
experiments Q. 

3 Conclusion 

We have carried out numerical simulations of a simplified model of magnetic grain in the 
thermal regime. We have carried two types of numerical simulations: switching time and 
switching field experiments. As far as switching time experiments are concerned, we found 
a gaussian behavior for the histogram P(lnr) of the logarithm of the relaxation times. This 
behavior is compatible with existing theories of slow relaxation phenomena. We find the 
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existence of a strong field regime, where both (In r) and the variance ((lnr- (lnr)) 2 ) tend to 
a constant as the temperature approaches zero. This regime was already predicted in ||, and 
called the strong field regime. If the magnetic field is smaller, the mean value and the variance 
diverge as the temperature approaches zero. It would be interesting to study the magnetic 
field h*, which separates the strong magnetic field region, as a function of the system size. 
It would also be interesting to carry out similar simulation for non Euclidian boundaries, 
in order to have a more precise idea of the importance of boundary effects. To do so, it is 
possible to study the case of a hyperbolic lattice, or a Cayley tree, with a finite fraction of the 
sites at the boundary even in the thermodynamic limit. Unfortunately, experimental data on 
magnetic grains are not available since these experiments require a large amount of data. We 
also calculated the streched exponent (3 as a function of T. Like in the case of experiments, 
we find that j3 may be superior to unity, and that (3 increases with T. 

We also considered the switching field experiment. There exists a master curve for the 
average of the switching field as a function of the switching rate and the temperature, which 
is consistent with the experimental data. 

The author acknowledges discussion with J.C Angles d'Auriac, B. Barbara, B. Dougot, K. 
Hasselbach, J. Souletie and W. Wernsdorffer. 
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Figure captions 



Figure 1: 

Snapshot of the spin configuration in a given numerical experiment, in the case of a 90x90 
square lattice. The magnetization is equal to 3/4. The inverse temperature is (5 — 1.6J, and 
the applied magnetic field is h = — J. A symbol denotes a down spin (parallel to the field) 
and the spin is up in the absence of symbol. The droplets are pinned at the corners. 



Figure 2: 

Distribution of P(lnr) as a function of lnr for different temperatures. The lattice is a disk 
with 2821 sites. The magnetic field is set to h = —1.5 from time t — 0. The different 
distribution and the fit to the Gaussian are plotted for T = 0.22J, 0.20J, 0.18J, 0.16J. 30000 
runs were carried out in each case. 



Figure 3: 

Mean value (lnr) as a function of temperature, for three values of the magnetic field h = 
—2J,h = — 1.5J and h = —J. The lattice is circular, with 2821 sites. 30000 runs were carried 
out in each case. 



Figure 4: 

Variance ((lnr — (lnr)) 2 ) as a function of temperature, for three values of the magnetic field 
h = —2J,h = —1.5 J and h = —J. The lattice is circular, with 2821 sites. 30000 runs were 
carried out in each case. 



Figure 5: 

Variance ((In r— (In r)) 2 ) as a function of the mean value (In r) , for three values of the magnetic 
field h = —2J,h = — 1.5J and h = —J. The lattice is circular, with 2821 sites. 30000 runs 
were carried out in each case. 



Figure 6: 

Variations of the streched exponent (3 with the temperature, and for h — — 1 and h = —1.5. 
The lattice is circular, with 2821 sites. 30000 runs were carried out in each case. 
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Figure 7: 

Master curve of the switching field experiment for the square lattice on the disk. We plotted 
{H sw )/Hg W as a function of f a (T,u). The Boltzmann constant is equal to unity, a = 1.8, 
Hg W = 2A,E = 5 and c = 1/25. On this curve, 7 series of measures at 7 different temperatures 
are superposed. The dashed line is the linear behaviour and corresponds to (H sw )/H® w = 
1 -/«(?>). 
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